A discriminating probe of gravity at cosmological scales 
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The standard cosmological model is based on general relativity and includes dark matter and 
dark energy. An important prediction of this model is a fixed relationship between the gravitational 
potentials responsible for gravitational lensing and the matter overdensity. Alternative theories of 
gravity often make different predictions for this relationship. We propose a set of measurements 
which can test the lensing/matter relationship, thereby distinguishing between dark energy/matter 
models and models in which gravity differs from general relativity. Planned optical, infrared and 
radio galaxy and lensing surveys will be able to measure Eg, an observational quantity whose 
expectation value is equal to the ratio of the Laplacian of the Newtonian potentials to the peculiar 
velocity divergence, to percent accuracy. We show that this will easily separate alternatives such as 
ACDM, DGP, TeVeS and f(R) gravity. 

PACS numbers: 



Introduction. — Predictions based on general relativ- 
ity plus the Standard Model of particle physics are at 
odds with a variety of independent astronomical obser- 
vations on galactic and cosmological scales. This failure 
has led to modifications in particle physics. By introduc- 
ing dark matter and dark energy, cosmologists have been 
able to account for a wide range of observations, from 
the overall expansion of the universe to the large scale 
structure of the early and late universe [ij . Alternatively, 
attempts have been made to modify general relativity at 
galactic 0] or cosmological scales A fundamental 

question then arises: Can the two sets of modifications 
be distinguished from one another? 

The answer is "No" if only the zero order expansion 
of the universe is considered. By allowing the dark 
energy equation of state wde to be a free function, 
the expansion history H(z) produced by any modified 
gravity can be mimicked exactly. Fortunately, struc- 
ture formation in modified gravities in general differs 
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general relativity. The difference we focus on here is 
the relationship between gravitational potentials respon- 
sible for gravitational lensing and the matter overden- 
sity. Lensing is sensitive to V 2 (0 — ip) along the line 
of sight where <f> and ip are the two potentials in the 
perturbed Friedman- Robertson- Walker metric: ds 2 = 
(1 + 2ip)dt 2 — a 2 (l + 20)g?x 2 and a is the scale factor. 
In standard general relativity (GR), in the absence of 
anisotropic stresses, <f> — —if), so lensing is sensitive to 



V 2 0. The Poisson equation algebraically relates V 2 </> to 
the fractional overdensity S, so lensing is essentially de- 
termined by 6 along the line of sight. This is a prediction 
of the standard, GR-based theory that is generally not 
obeyed by alternate theories of gravity. 

Testing this prediction is non-trivial. Astronomers of- 
ten use the galaxy overdensity as a probe of the underly- 
ing matter overdensity, but the two are not exactly equal. 
Here we propose a test of this prediction which is rela- 
tively insensitive to the problem of galaxy bias. The basic 
idea is simple: 

• Extract the matter overdensity at a given redshift 
by measuring the velocity field. Matter conser- 
vation relates velocities to the overdensities. The 
measurement of the velocity field can be accom- 
plished by studying the anisotropy of the galaxy 
power spectrum in redshift space. 

• Extract the lensing signal at this redshift by cross- 
correlating these galaxies and lensing maps recon- 
structed from background galaxies. 

More quantitatively, the galaxy-velocity cross power 
spectrum P g g = — {5 a (k)9(— k)} can be inferred from 
redshift distortions in a galaxy distribution. Here, = 
V • v/H(z) and v is the comoving peculiar velocity. In 
the linear regime, matter conservation relates to 5 by 
= —8/H = —(35, where (3 = d\nD/d\na and D is 
the linear density growth factor. So, P g g = f3P g s, sat- 
isfying the first goal above. Cross correlating the same 
galaxies with lensing maps constructed from galaxies at 
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higher redshifts, P^2 



can be measured. The ra- 



tio of these two cross-spectra therefore is a direct probe 
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of V 2 (</> — ip)/(/38). It does not depend on galaxy bias 
or on the initial matter fluctuations, at least in the lin- 
ear regime. Modifications in gravity will in general leave 
signatures in either (3 and/or the Poisson equation. 

Galaxy- Velocity Cross-correlation. — A galaxy's 
peculiar motion shifts its apparent radial position from 
x z to x z = x z +v z / H(z) in redshift space, where v z is the 
comoving radial peculiar velocity. The coherent velocity 
component changes the galaxy number overdensity from 
S g to 5g ~ 5 g — V z v z j H(z). Galaxy random motions mix 
different scales and damps the power spectrum on small 
scales. The redshift space galaxy power spectrum there- 
fore has the general form ([17J and references therein) 



P g S (k) = [Pg(k) + 2u 2 P g g(k) + U 4 Pg(k)] F 



H 2 (z) 



(1) 



where u = kn/k is the cosine of the angle of the k vector 
with respect to radial direction; P g , P g g , Pg are the real 
space galaxy power spectra of galaxies, galaxy-6* and 9, 
respectively; a v is the ID velocity dispersion; and F(x) 
is a smoothing function, normalized to unity at x = 0, 
determined by the velocity probability distribution. This 
simple formula has passed tests in simulations on scales 
where 8 < 1 The derivation of Eq. JTJ) is quite 

general, so it should be applicable even when gravity is 
modified. 

The distinctive dependence of Pg on u allows for si- 
multaneous determination of P g , P g g and Pg 18]. The 
parameters we want to determine are the band powers of 
P g g(k) 1 defined such that P(k) = P a if k a < k < k a+ \, 
where ki < k 2 < ■ ■ ■ < k a < ■ ■ ■ . We denote Pa^ 
as the band power of P g g. For a in each k bin, 
we have a measurement of Pg, which we denote as Pj. 

The unbiased minimum variance estimator of Pa is 



1 Distance D and H are required to translate the observed galaxy 
angular and redshift separation to k. In general, errors in D and 
H measurements cause both horizontal and vertical shifts in the 
Eq plot. Both D and H will be measured by methods like type 
la supernovae and baryon acoustic oscillations with 1% accuracy, 
much smaller than the k bin size adopted, so the horoziontal 
shift is negligible. Errors in D show up in both P g g and the 
C Kff — » ^v 2 (i/>— 0)9 inversion through I = kD and thus largely 
cancel in evaluating Eq. Errors in H(z) only show up in P g g 
measurement and thus cause a net shift in the value of Eq. For 
1% error in H, the fractional error in Eq is < (n c g + 3)1% < 
3%. Here, n e ff is the effective power index of the corresponding 
power spectra. For the fiducial ACDM cosmology, it is negative 
in relevant k range. Thus errors induced by uncertainties in 
D(z) and H(z) measurement will be sub-dominant, except for 
SKA, which requires better control over systematic errors in D 
and H measurement. For simplicity, we neglect this potential 
error source. Measuring P g g also requires to marginalize over 
<t„. However, in the linear regime k < 0.2fe/Mpc, k 2 al/H 2 < 1 
and F{k' 2 u 2 ul/H 2 ) ~ 1, for typical value u v ~ 300 km/s. Thus 
the exact value of cr v is not required for our analysis. Without 
loss of generality, we adopt cr v = 300 km/s. 



P = E WiPt, where W, = ^r(Ai + A 2 u? + X 3 uf). Here, 

Fi = F(kui<r v /H), of is the variance of Pi and the three 
Lagrange multipliers A Q (a = 1,2,3) is determined by 



A = (0,i, 0)-A- x 
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Galaxy-galaxy lensing. — Weak lensing is sensitive 
to the convergence k, the projected gravitational poten- 
tial along the line of sight: 



V 2 (iP-0)W( X ,Xs)d X 



(3) 



Here, W is the lensing kernel. For a flat universe, %, \s 
are the comoving angular diameter distance to the lens 
and source, respectively. Eq. [3] is a pure geometric result 
and can be applied to any modified gravity models where 
photons follow null geodesies. 

A standard method to recover the lens redshift infor- 
mation is by the lensing-galaxy cross correlation. For 
galaxies in the redshift range [21,22], the resulting cross 
correlation power spectrum under the Limber's approxi- 
mation is 



C Kg (l) = 4 / n g (x)d X 
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Here, Xi,2 are the comoving angular diameter distance to 
redshift 2i )2 and x 1S the mean distance. The band power 
Pi 2] of P v 

2 (4>-<t>) g is defined at the same k range as P^\ 
In practice, we measure the band power C Kg (l, A/), cen- 
tered at I with band width Al. The weighting f a (l, Al) 
is defined correspondingly. For each I, only a fraction of 
a having f a (l, Al) ^ contribute. 

A discriminating probe of gravity. — With the 
above measurements, one can construct an estimator 



E G = 



C Kg (l,Al) 



3H$a-ij: a f a (l,Al)P< 



whose expectation value is 

(Eg) 
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The fractional error on Eg is 

AC 2 , Y, a Il((5P^Y) 
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LAM0ST/AS2+LSST 




ADEPT+LSST , 3<z<1 7 
SKA 



SKA 0.85<z<1.15 



SKA 1.8<z<2.2 
, I 



TABLE I: Summary of target surveys. 
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FIG. 1: Eq as a smoking gun of gravity. Error estimation is 
based on ACDM and error bars are centered on the ACDM 
prediction (black solid straight line). We only show those k 
modes well in the linear regime. For clarity, we shift the error 
bars of LAMOST/AS2+LSST and ADEPT+LSST slightly 
rightward. Irregularities in the error-bars are caused by irreg- 
ularities in the available discrete k modes of redshift distor- 
tion. Dotted lines are the results of a flat DGP model with 
tto = 0.2. Dashed lines are for f(R) = -Ai# 2 exp(-_R/A 2 -ff 2 ) 
with A2 = 100. Differences in expansion histories of these 
models are of percent level at 2 < 2 and are not the main 
cause of differences in Eg- Solid lines with wiggles are for 
TeVeS with Kb = 0.08,0.09,0.1, where the lines with most 
significant wiggles have Kb = 0.1. 



where AC 2 = [C 2 Kg + (C K + C»)(C g + C» )]/(2/A// sky ). 
Here, C K , , C gi C g are the power spectra of weak 
lensing convergence, weak lensing shot noise, galaxy and. 
galaxy shot noise, respectively, and / s ky is the fractional 
sky coverage. Errors on Eq at any two adjacent bins are 
correlated, since they always share some same k modes. 
However, by requiring l a /xi = la+i/x^, where h < h < 
■ ■ ■ < l a < ■ ■ ■ and k a — I a /X2, Eq measurement at each 
I bin only involves two k bins and thus only errors in 
adjacent bins are correlated. 

We choose ongoing/proposed spectroscopic surveys 
LAMOST, AS2, ADEPT and SKA as targets of redshift 
distortion measurements, and LSST and SKA as targets 
of lensing map reconstruction. SKA lensing maps can 
be constructed through cosmic magnification utilizing its 
unique flux dependence, with S/N comparable to that of 
LSST through cosmic shear [191 ] . Survey specifications 
are summarized in TABLE I. The fiducial cosmology 
adopted is the ACDM cosmology, with the WMAP best 
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fit parameters Qq = 0.26,17a = 1 — Ho,h = 0.72, erg = 
0.77 and n s = 1. The result is shown in figure [TJ In gen- 
eral, at k < 0.1/i/Mpc, cosmic variance in C Kg and P g g 
measurements dominates the Eq error budget, result- 
ing in decreasing error-bars toward larger k. This makes 
/sky and the lensing source redshifts the two most rele- 
vant survey parameters for Eq error estimation. Since 
systematic errors in LSST photometric redshifts can be 
controlled to better than 1%, errors in Eq measurements 
of LAMOST/ AS2+LSST and ADEPT+LSST caused by 
source redshift uncertainties are sub-dominant. 

We restrict our discussion to sub-horizon scale pertur- 
bations and express equations hereafter in the Fourier 
form. Four independent linear equations are required 
to solve for four perturbation variables 5, 9, ip and <fi. 
The mass-energy conservation provides two: S + H 9 = 
and HO + 2H 2 9 - k 2 ^/a 2 = 0. For at least ACDM, 
quintessence-CDM, DGP and f(R) gravity, the other two 
takes the general form 



(j) = -r)(k,a)ip , 



(8) 



Here flo is the cosmological matter density in unit of the 
critical density p c = 3H^/8nG. Refer to 0, [U, for 
other ways of parameterizations. MOND has extra scalar 
and vector perturbations and does not follow the general 
form of Eq. 0,0. 

(1) ACDM: 77 = 1, G eff = 1 and E G = Q /f3. Dynam- 
ical dark energy will have large-scale fluctuations [-ioj ]. 
Furthermore, it may also have non-negligible anisotropic 
stress and is thus able to mimic modifications in gravity 
2l|. But, for models with large sound speed and neg- 
ligible anisotropic stress, such as quintessence, these are 
negligible at sub-horizon scales and Eq. 03 still holds. 

(2) Flat DGP: V = [1 - l/3/3 DG p]/[l + 1/3/W], 
G off = 1 and E G = Q Q //3, where /3 DG p = l-2r c fT(l + 
H/3H 2 ) < and r c = H /(l - Q ). Q a differs from that 
of ACDM, in order to mimic H{z) of ACDM. 

(3) f(R) gravity: in the sub-horizon limit, G c g — 
(1 + IrY 1 Q and 77 = 1 [HI, with f R = df/dR\ B where 
B denotes the FRW background. This falls naturally out 
of a conformal transformation of the expression for Eq 
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in the Einstein frame into the Jordan frame, noting that 
Einstein frame scalar field fluctuations are negligible on 
sub- horizon scales [l2} • We numerically solve the full per- 
turbation equations in the Einstein frame since it is com- 
putationally simpler [12$ and then conformally transform 
to the Jordan frame, which we choose as the physical 
frame, evaluating (3 such that Eq = Qo/(l + /r)/3. In the 
limit that f R -> 0, e.g. for f(R) ~ Ai£T 2 exp(-i?/A 2 i^) 
111 ] with Ai «C A2, the evolution is observationally equiv- 
alent to ACDM. For modes that entered the horizon prior 
to matter-radiation equality, as we consider here, /3, and 
therefore Eq, is scale invariant for IR modifications to 
gravity, with fn > 0. 2 The scale independence of Eq 
holds in ACDM, Quintessence- CDM and DGP. An ob- 
served scale-independent deviation in Eq from ACDM 
could signify a special class of modified gravity, as shown 
in Fig. Q] 

(4) TeVeS/MOND. Besides the gravitational met- 
ric, TeVeS Q contains a scalar and a vector field. These 
new fields act as sources for the gravitational potential 
<t> in the modified Poisson equation and can change the 
evolution of cosmological perturbations with respect to 
standard gravity 0, 0] . We considered a TeVeS model 
with fij = 0.05, VL V = O.f 7, £1a = 0.78 and we adopted a 
choice of the TeVeS parameters that produces a signifi- 
cant enhancement of the growth factor. The TeVeS Eq 
is significantly different from the standard Eq (Fig. [I}. 3 
It exhibits scale dependence with accompanying baryon 
acoustic wiggles. Both features are due to the vector 
field fluctuations, which play a significant role in struc- 
ture formation [7J. These fluctuations decrease toward 
small scales and cause the scale dependency of Eq. We 
also checked that they affect the final shape of the acous- 
tic oscillations of the other components significantly. As 
a result, oscillations in </>, ip and 5 do not cancel out per- 
fectly in TeVeS when we take the ratio, thus producing 
the wiggles in Eq. 

For the four gravity models investigated, differences 
in Eq are much larger than observational statistical un- 
certainties. Planned surveys are promising to detect per- 
cent level deviation from GR and should distinguish these 
modified gravity models unambiguously. 



2 Scales larger than the horizon at matter-radiation equality are 
suppressed [T3l and, if measurable, would have a scale dependent 
increase in the value of Eq in comparison to the small scale value. 

3 To simplify the numerical treatment of the TeVeS perturbations 
equations while retaining a good qualitative description of all the 
significant physical effects at the same time, we introduced sev- 
eral approximations. Namely we assumed instantaneous recom- 
bination and employed the tight coupling approximation between 
baryons and photons at all scales before decoupling; moreover 
we evolved perturbations in the massive neutrino component in 
a simplified way by switching off neutrinos perturbations when 
they were below the free steaming scale and treating them as a 
fluid above the free streaming scale. 



At large scales, gravity is the only force determining 
the acceleration of galaxies and dark matter particles. So 
we assumed no galaxy velocity bias. As statistical errors 
in Eq measurements reach the 1% level (Fig. I), sev- 
eral systematics, besides the one discussed in footnote 1, 
may become non- negligible. One is the accuracy of the 
redshift distortion formula (Eq. [T]) , which may be prob- 
lematic for those modes with large u, even at very linear 
scales 17j • A remedy is to exclude them when extracting 
P g e , at the expense of statistical accuracy. As discussed 
before, accuracy of Eq measurement is dominated by 
accuracy of fv 2 (v>-</>)ff measurements and is thus less af- 
fected. A less severe one is the nonlinear evolution, which 
becomes non-negligible where the matter power spectrum 
variance ^0.1. In general relativity, nonlinear cor- 
rections to density and velocity differ (Fig. 12, |22[). A 
direct consequence is that Eq develops a dependence on 
the matter power spectrum. Similar effects in modified 
gravity models are expected. This can be corrected by 
high order perturbation calculations, which should work 
well where < 0.2. 
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